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A real-space-renormalization method for the frequency dependent conductiv- 
ity of the periodic approximants of the Fibonacci chain is developed. This 
scheme is based on the known 2x2 transfer matrices and additional 5x5 
matrices which allow an efficient numerical evaluation of the Kuho formula. 
Numerical results are presented. 
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1. INTRODUCTION 

Quasicrystals (QC)0 and their periodic approximants reveal unusual 
transport properties, e.g., extremely low conductivity for alloys of metallic 
constituenta3, and it is generally believed that these properties are directly 
related to the long range but nonperiodic order. However, no coherent the- 
oretical framework has been successfully established. For instance, previous 
calculations of the density of states (DOS) revealed a set of peaks and gaps 
of small width (10 — 100meV)i but recent studies nrovide convincing argu- 
ments that these structures are numerical artifacta3. Most of the previous 
numerical studies of the conductivity were based on the Landauer formula p 
but only few studies were made in the framework of linear response theoryt . 

One-dimensional QCs, such as the Fibonacci chain (FC), have been 
studied in great detail by many researchers, because exact results and an- 
alytical treatments are possible. The most striking feature of this model 
is that quasiperiodicity induces long range correlations giving rise to an 
intermediate state of localisation. These so-called 'critical states' are associ- 
ated to singular continuous spectra of the Hamiltoniani, while all states are 



*To appear in Journal of Low Temperature Physics, March 2002 



2 



D. Walther and R. v. Baltz 



ep ea e-y ep ep e-, ep 

th ts tL tL ts tL ts tL tL ts tL th ts 

Fig. 1. Fifth generation standard Fibonacci lattice = LSLLSLSLL 
SLLS (/s = 13) corresponding to the FC supercell of a periodic approxi- 
mant. 



locahzed for randomly disordered systems in 1-D. Most progress in investi- 
gating 1-D models has been made by methods based upon the trace-map 
and real-space decimation techniques. To the best of our knowledge nothing 
comparable has been worked out for the current-current correlation func- 
tion which allows for a renormalization group approach of the frequency 
dependent conductivity. 

Our paper is organized as follows. Sec. 2 summarizes the physical model, 
in particular, the relation between the geometrical and the electronic struc- 
ture which has been worked out in a preceding paperQ. In Sec. 3, a suitable 
expression of the Kubo conductivity is derived in terms of the known 2x2 
transfer matrices and additional 5x5 matrices A, which allow to set up a 
real-space renormalization scheme and its numerical implementation. Sec. 
4 gives some numerical results and conclusions. 



2. FIBONACCI CHAIN 

The FC is specified by two different types of "bonds" denoted by L 
("long") and S ("short"), see Fig. 1. One may describe the geometric struc- 
ture of such lattices by words w{L, S), i.e. strings in the symbols L, 5 rep- 
resenting the corresponding linear arrangement of the "letters" L and S. 
These words can be generated by an (infinite) repetition of the recursion 
law (Nielsen transformation (NT)) 

Wn = Wn-l * 'Wn-2, W-i = S, Wq = L , (1) 

where the symbol * is defined as the concatenation of two strings and n = 
1, 2, . . . The length of the FC in the n^^ generation \wn\ = fn, i-e. the sum of 
symbols L and S in Wn{L, S) satisfies the recursion relation for the Fibonacci 



fn — fn-l + fn-2, /-I — /o — 1 • 



(2) 
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Alternatively, the n order approximants can be generated by substitu- 
tions (morphisms), which operate on the symbols L, S rather than on words 
WniL, S) 



The mapping C may be viewed as the orbit of L, i.e. Wn = C^{L): 

L^LS^ LSL ^ LSLLS ^ . . . ^ LSLLSLSLLSLLS ■ ■ ■ 
In our study, we employ the following tight-binding Hamiltonian 



where \fi > are Wannier states centered at sites fi & G = [—M,N — 1] C 
Z, and the real e^,t^^jy denote the site-energies and the nearest-neighbour 
transfer integrals, respectively, {t^^i, > 0). The FC has the speciality that 
there is no S-S bonding, hence only three of the possible four combinations 
occur: 



To each word w = {y^)^,eG = V^mV-m+i ■■■Vn^ with £ {L, S}, M + 
N > 2, we uniquely assign a dual word S = {a^)^^\j = {o-_j^j, . . . , Cjy-i)! 
where each symbol € {a, (3,^} will be related to the pair by 
the map (LL, LS, SL) <-> (a,/3, 7). For example, the dual word pertaining 
to W5 is S'^ = 'yP'jaP'yP'-faP'ya, with G = [0, 12], see Fig. |l|. Let be = 
(crp^zi^ and S^, = ('7p^=^p the dual words of u, and v respectively. We 
define the product of two dual words, corresponding to the product w = u*v, 
by A =: T,^, with = (cJ^^, • • • ctJ^_i, ct^, cr!!^, ■ • • ct^.J, ct„ := 
{y^,y'^p), where yj^ (y-p) is the last (first) letter of u {v). 

To set-up a renormalization scheme we approximate the aperiodic FC 
by a periodic approximant of length 2N (M = N), 



= S„AS„A... ASn^ , G=[-N,N-1], N = N'f„, (6) 



where /„ — 1 denotes the number of symbols o"^ of S^. = ^ , t^+Zn = , 

(J, = fi + mf^ where jl € [0, /,j — 1] labels the 'atoms' inside the Fibonacci 
supercell and m G [—N', N' — 1] labels the vectors of the Bravais lattice 




(3) 



?i = ^ l/Li > < /ij + ^ 1^ > t^^i, < V 



(4) 




(5) 
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of the periodic approximant. Obviously, there is a unique correspondence 
between the dual word S and the FC Hamiltonian Ti^ . 

In site representation, the eigenstates of Eq. @ can be represented 
in terms of the two fundamental solutions {P^} and {P^} of the discrete 
Schrodinger equation 

{n^y)^ := + t^^^^-iiE) + e,^,iE) = E^^{E) , (7) 

such that the (unnormalized) sokitions are 

^iiE):=V,{E)+mi{z)P,{E). (8) 

Pfj_{E), Vfj,{E) respectively denote polynomials of degree ^ and /i — 1 with 
initial conditions P_-^ = 0, Pg = 1/tg and 'P_i = — l/tg, "P,, = iS. In addition, 
we have to fulfill the Bloch conditions ^^±/„ = exp(ibzA;a)^'^, where a = /„ 
is the "lattice constant" and k = k{E) is the wave vector. According to 
time-reversal symmetry a state of energy E' is a twofold degenerate which 
will be denoted by \E,±) and ^(-E) = (a*|±), respectively. m±{E) denotes 
the Titchmarsh-Weyl functionll3 

i^{E)+t,V^^_AE) 
m±iE) = — , (9) 




Ci{E) := e^^'^^>^- =^±^^-fj -1, (10) 

\ = *o(^.„-^,„-J, (11) 

where f„k{E) = {pt^)/N' = arccos(A^^ {E)/2) G [0,7r), and < p < N' . 

For the standard one-atomic tight binding chain with equal t' s and 
e' s A^. {E) and P^{E) = V^-i are, apart from a scale factor , Chebyshev 
polynomials T^{x), U^{x) of first and second kind, respectively^. ^^{E) = 
g±*Mfc^ ^(A:) =e-2tcos(A:). 



3. CONDUCTIVIY AND REAL SPACE 
RENORMALIZATION 

In the linear response theory the real part of the conductivity is given 

byS 

»e{.<'.)M} = r m^M±Mr<-^^E..)^E. (12) 

r(")(E,cj) = Tr{J6{E -'H^^„))JS{E + nw -n^in^)} , (13) 
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Fig. 2. Density of states (left), bandstructure (middle), and (unrestricted) 
solutions of Eq. (^) (right) as function of frequency. 6*'* order generation 
FC with ea = e/3 = = and = 1.1, t^ = 1.0. 



where f{E) is the Fermi function and J is the current operator 
i ieh 



m 



j^lu\^J■){J^\, jfj.u = -t^,+id^u-i+tf^5^iu+i- (14) 



At zero temperature and after some nontrivial algebra Eq. ([l2| ) becomes 



g2 ?Tr|r0T0 + T ^0T"^0 - 200| 



9,,, I A' - A' I . /a _ (' A ('FS\^ 



(15) 



E = 
i^>0 



with A' := dA/d^; and ? = sign(P^^_j(£;)Pj.^_j(^ + ^o;)), respectively. 

(Note, in Eq. ([l5|) we have omitted the dc-part g^^^6{uj) which vanishes 
in the limit n ^ oo). All quantities are functions of E and/or E + f\w and 
@{E,E + huj) = @{E + huj,E),T{E) = T{E + huj), etc ... . 

-Eoj denotes the energy of an initial occupied (lower) band from which 
transitions to an empty (upper) band can take place. By conservation of 
energy and (crystal-) momentum E = E^ is a solution of 

k{E + huj) = k{E) A^, (E + huj)- A^ (^) = . (16) 



Remarkably, the set of all solutions of Eq. (p^q), denoted by E{uj), also in- 
cludes unphysical solutions which smoothly join E^, when both E + huj and 
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E are in the gap regions (where k is imaginary), see Fig. |2|. Of course, only 
the subset in the spectrum of "H^, should be included in Eq. (]T5|). 
To evaluate Eq. (|l5|) we introduce a matrix notation 



(17) 



P/,"^ are the polynomials associated to which are obtained by replacing e^, 

by Cf^+u, t^+u, e.g. = P^'Ji . In addition, we define @{E, E+huj) = 
with 



E 



2U 



12 



A' 



(n) 



^(n) 



(18) 



A^"^ are components of the 5 x 5-matrix 



d"'\E,E + hw) 


= A/„-2 • Af„ 


-3 • . . . • 


Ai, 




(19) 


where 
















fSf,^^, + 1 














1 
















A^(E,^ + M = 





1 






































-(Am 


















(20) 


We have set = {E 




,Sfj^ = {E + huj - 








and 






Then = {C^s 









(2C^_it^) and 



Qn'Tfj.-i/ In contrast to the matrix 

1 



"f -1 



■T^^, withT,^ 







1 

E-e 



(211 



A(T;j_i(T;jo-;j+i ) see 



which corresponds to each site o"^ (i.e. to each "dual letter", such that = 
TfjJ) uniquely, the correlates three sites (i.e. a triple of three dual letters) 
C7^_i, cr^ and cr^+i, respectively. Thus we may write A^ 
Fig. |. 

The self similarity of EC shows up in the structure of Eqs. (1£,21) which 
will lead to a real-space renormalization scheme of the structure of the trace- 
map methodlli or path renormalization used for the calculation of the local 
DOS@. To set up the renormalization scheme we consider the NT of words, 
cf. Eq. d). 



Xl-s,s{Wk-l,Wk-2) = {Wk,Wk-l) , Wk = {Wk-lY Wk-2 [Wk-l] 



(22) 



Conductivity of the Fibonacci Chain 



7 









i 












7 


/3 7 




fi 7 fi 7 a 




7 




^^ 7 


/3 




Tg * T. t 


• 

To 


• • • • • - 


* Ts 


* T. * 


• 

Ta * 






; 






i 












7 


13 7 


a 


13 1 13 ~i a 


13 


7 




P 7 








-•- 


• • • • • - 


Q^*Ao 


;37*A,37Q 


• — 


*Aq^t 










AMI 


55 




AI3) 







Fig. 3. Example of the RTs (24) and (Bn) for a 5 generation FC. Arrows 



indicate the dual letters a = a)"'' as given by Eq.p5|). Ti * T2 * T; 

Jn 

T3 • T2 • Ti. etc., where the dot denotes the usual matrix-product. 



with s = 0.1, and the initial words w-i = S, wq = L. Next, we write the 
trace-mapEj in the form (recall, that A^. = Tr{T}) 



where L = Xi^, R = Xq^i. Note, is the renormalization transforma- 
tion (RT) corresponding to the NT L etc.. Furthermore, the quantity A = 
Af +Af +Af —A/, ^Af, , A/, is invariant under the transformation 

so that the trace-map Eq. ( p3[ ) represents a dynamical system which is 
essentially two-dimensional. Therefore, (A/^^^ , A/^ ) = (^i^)"(A/2 , A/^ ). 
In addition, the NTs (^) implies 



(24) 



By inspection of Eq. (21) we recognize that T^"^ does not correspond to 
S„. but to the cyclic permutated word T(S„) = {a-^, . . . ,a, ), a. = cr^. 
The corresponding word idlll w^^\L,S) = TZ o (C)^~^{L) which implies 
(r("+i),T(")) = («L)""^o«i?(T(2),T«). 

The case of the matrices A^"-* is a little bit more involved, see Fig. |^. 
Let A("-i) and A^"^^) correspond to the dual words and 
respectively. Then the NT L implies i(S„_i,S„_2) = i'^n,'^n~i), where 



^ J n — l •> •> 71 — 2 



(25) 
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Eq. dH) corresponds to A^") = A^"-^) ^ ^(n-i) — A^""^) .^^.^(n-i)^ ^^^^^ 
we have introduced a "sandwich" -matrix 



5n = A, 



1) • A 



(n-2) (n-2) 
Jn-2 Jn-2 ^ " 



As a result, the NT L leads to the following RT for the A^'') 

^ = / '^T'^T ■ '^'37/3 , even 
^ \ Aa/3^ • A^a/3 , k odd ' 



(26) 

(27) 
(28) 



so that (a("+^\ A(")) = (^I.)'^(A(2), a(^)) holds. Eqs. (|1,|1|1) provides 
a simple renormalization scheme to calculate the frequency dependent con- 
ductivity of the standard FC. This is our main result. 



4. RESULTS AND CONCLUSIONS 

Figs. |2|, 1^, ^ display some numerical results for the DOS, minHaandstruc- 
ture, transition energies, and the conductivity. We note, that the expression 
of a^'^^Lv) in Eq. ( |l5|) is nonzero only for ui > 0. For w — > 0, merely a^^\0) 
contributes because none of the lines {eI^^ , u) cut the axis a; = for energies 
within the spectrum. Therefore, a^"'\uj) vanishes in the limit w — > 0. (In 



addition, the dc conductivity (see after Eq. (15)) vanishes in the limit of 
infinite aperiodic FC, n — > oo.) This is characteristic for the (infinite) TB 
model, as H^^^ is a bounded operator. Thus, there is always a gap in the 
neighbourhood of a; = 0. If one introduces some mechanism of dissipation 
or disorder one may expect that the small frequency-gap revealed in Fig. ^ 
may be smeared out (leading to metallic-like behaviour), but will persist in 
the case of Fig. ^ (semiconductor-like behaviour). 

The main numerical challenge turned out to be the calculation of the 
energies S^, i.e., the zeros of Eq. ([l6|). This task implicitly implies polyno- 
mials of very high order. Such operations are numerically unstable biiLcould 
be managed with the help of Bailey's multiprecision software packetli3. 

Starting from the self-similarity of the FC we found a real space renor- 
malization scheme for the Kubo-conductivity which is suitable for numerical 
implementation. Hereto, this scheme stands alongside with the powerful tool 
of the trace map for the spectrumlll and the path renormalisation scheme for 
the local DOSQ. For details and generalizations of our approach to arbitrary 
morphisms (like Thue-Morse, Rudin-Shapiro, period doubling) and n-letter 



alphabets, we refer to Ref. 16 
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Fig. 4. Density of states N{E) (left) and conductivity (t{uj) (right) of the FC 
with a supercell of /n = 233 sites. Fermi energy ep = 0, ts = 1.0, = 1.5, 
and = 6/3 = = (arbitrary units). 




Fig. 5. Density of states (left) and conductivity (right) of the same FC as in 
Fig. 1^ but with different parameters: ep = —0.2, — = e^j = e-^ = 1.0. 



ACKNOWLEDGMENTS 

This work was supported by the Deutsche Forschungsgemeinschaft DFG 
through the "Schwerpunktsprogramm" Quasikristalle. D. H. Bailey and 
NASA Ames are kindly acknowledged for making available their Multiple 
Precision Floating Point Computation Package^. 



10 



D. Walther and R. v. Baltz 



REFERENCES 

1. For an introduction see, e.g., C. Janot: Quasicrystals: A Primer, Clarendon 
Press (1994); F. Axel and D. Gratias, eds.. Beyond Quasicrystals, Springer 
(1995). 

2. E.g., F. Hippert und D. Gratias, eds.. Lectures on Quasicrystals, Les Editions 
de Physique, Les Ulis (France) (1994); Z. M. Stadnik, ed.. Physical Properties 
of Quasicrystals, Springer Series in Solid State Sciences 126, Springer (1999). 

3. T. Fujiwara, Phys. Rev. B 40, 942 (1989); J. Hafner and M. Krajct, Phys. Rev. 
B 47, 11 (1993); T. Fujiwara, S. Yamamoto, and G. T. de Laissardiere, Phys. 
Rev. Lett. 71, 4166 (1993); D. Mayou, C. Berger, F. Cyrot-Lackmann, T. Klein, 
and P. Lanco, Phys. Rev. Lett. 70, 3915 (1993). 

4. E. S. Zijlstra and T. Janssen, Phys. Rev. B 61, 3377 (2000); E. S. Zijlstra and 
T. Janssen, Europhys. Lett. 52, 578 (2000). 

5. B. lochum and D. Testard, J. Stat. Phys. 65, 715 (1991); W. Salejda and P. 
Szyszuk, Physica A 252, 547 (1998); B. Sutherland and M. Kohmoto, Phys. 
Rev. B 36, 5877 (1987). 

6. S. I. Ben-Abraham and A. Joseph, Proceedings of the 5th Int. Conf. on Qua- 
sicrystals, C. Janot and R. Mosseri, eds.. World Scientific (1995), p. 621; A. 
Lahiri, Phys. Rev. B 53, 3702 (1996); S. Roche and D. Mayou, Phys. Rev. Lett. 
79, 2518 (1997). 

7. D. Walther and R. v. Baltz, Phys. Rev. B 55, 8852 (1997). 

8. Leonardo Pisano (1170-1250), (named Fibonacci, from filius Bonacci), 
Liber abbaci (1202). For a historical review see: http: / /www- groups. dcs.st-| 



and .ac.uk/ "-history /Mat hematicians /Fibonacci . html 



9. E.g., C. Brezinski, Pade-Type Approximation and General Orthogonal Poly- 
nominals., Birkhauser-Verlag (1980). 

10. H. Weyl, Ann. Math. 68, 220 (1910); E. C. Titchmarsh: "Eigenf unction 
Expansions Associated with Second-Order Differential Equations^\ Oxford, 
University Press (1946). 

11. Natl. Bur. Stand. Appl. Math. Ser. No 55, edited by M. Abramowitz and I. A. 
Stcgun, (U.S. GPO, Washington, DC) Handbook of Mathematical Functions, 
Dover (1968). 

12. E.g., A. Douglas Stone und Aaron Szafer, IBM J. Dev. 32, 384 (1988); H. U. 
Baranger und A. D. Stone, Phys. Rev. B40, 8169 (1989). 

13. M. Kohmoto, L.P. Kadanoff, and C. Tang, Phys. Rev. Lett. 50, 1870 (1983); 
S. Ostlund, R. Pandit, D. Rand, H. J. Schellnhuber, and E. Siggia, Phys. Rev. 
Lett. 50, 1873 (1983). 

14. See the Nielsen-Graph, Fig. 5, of Ref. |^. Note, the sequence of substitutions 
C, IZ is reversed with respect to the NTs L, R. Here, the path which connects 

and a[''^ = (3['''' provides the sequence of NTs of the RG. 

15. D. H. Bailey, Multiprecision Translation and Execution of Fortran Programs, 
ACM Transactions on Mathematical Software, 19, no. 3, Sept. 1993, p. 288-319. 

16. D. Walther, Dissertation, Universitat Karlsruhe, Nov. 2001. 



